Reconstruction of the deceleration parameter and the equation of state of dark energy 
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The new 182 gold supernova la data, the baryon acoustic oscillation measurement and the shift 
parameter determined from the Sloan Digital Sky Survey and the three-year Wilkinson Microwave 
Anisotropy Probe data are combined to reconstruct the dark energy equation of state parameter 
w(z) and the deceleration parameter q(z). We find that the strongest evidence of acceleration 
happens around the redshift z ~ 0.2 and the stringent constraints on w(z) lie in the redshift range 
z ~ 0.2 — 0.5. At the sweet spot, —1.2 < w(z) < —0.6 for the dark energy parametrization w(z) = 
u>o + w a z/(l + z) 2 at the 3a confidence level. The transition redshift Zt when the Universe underwent 
the transition from deceleration to acceleration is derived to be z t — 0.36^qq|. The combined data 
is also applied to find out the geometry of the Universe, and we find that at the 3a confidence level, 
|fifc| J$ 0.05 for the simple one parameter dark energy model, and —0.064 < fit < 0.028 for the 
ACDM model. 

PACS numbers: 98.80.-k,98.80.Es 
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I. INTRODUCTION 

The discovery of the accelerated expansion of the Uni- 
verse by the supernova (SN) la observations [H imposes 
a big challenge and provides opportunities to theoreti- 
cal physics. The more accurate SN la data @, 0, H, 
the three-year Wilkinson Microwave Anisotropy Probe 
(WMAP3) data [|, and the Sloan Digital Sky Survey 
(SDSS) data [| tell us that the Universe is almost spa- 
tially flat, and dark energy (DE) with negative pressure 
contributes about 72% of the matter content of the Uni- 
verse. Although the existence of DE were verified by 
different observations, the nature of DE is still a mystery 
to us. For a review of DE models, one may refer to Ref. 
0- 

Many parametric and non-parametric model- 
independent methods were proposed to study the 
evolutions of the deceleration parameter q{z), the 
DE density £Ide(z)> the DE equation of state 
(EoS) w(z} x and the geom etry of the Universe 





In the reconstruction of q(z), it was found that the 
strongest evidence of acceleration happens at redshift 
z ~ 0.2 [1, 0, [l(|, and the evidence of the current 
acceleration is not very strong [9( and model dependent 
[To| . Previous studies on the reconstruction of w(z) 
also showed that the stringent constraint on w(z), or 
the swe et spot, happened around redshift z ~ 0.2 — 0.5 
d El El El EE EH- As Riess et al. pointed out, 
the use of additional parameters to reconstruct w{z) 
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does not provide a statistically significant improvement 
on the fit of the redshift-magnitude relation, so we 
discuss one- and two-parameter models only. In this 
work, we first use the simple two-parameter model 
q(z) = 1/2 + (qiz + q 2 )/(l + z) 2 [10j to reconstruct q(z), 
then we use the three popular two-parameter models 
w(z) = w + w a z/ (l + z) w{z) = w + w a z/ (1 + z) 2 
IB and Q DE = l-n m -A 1 -A 2 + A 1 (l + z) + A 2 {l + z) 2 
14j to reconstruct w(z). The purpose of this work is 
to see if the stringent constraints on q(z) and w(z) still 
happen around z ~ 0.2 — 0.5 when we use the new 
182 gold SN la data compiled in Q|. The geometry 
of the Universe is also discussed by fitting the simple 
one-parameter model w(z) = w exp[z/(l + z)]/{\ + z) 
El to the combined SN la, SDSS and WMAP3 data. 



This paper is organized as follows. In section II, we 
study the property of q(z) by fitting the parametriza- 
tion q(z) = 1/2 + (qiz + <z 2 )/(l + z) 2 to the new 182 
gold SN la data compiled in [4(. In section III, we apply 
the popular parameterizations w(z) — wq + w a z/(l + z), 
w(z) — wq + w a z/(l + z) 2 and flnE = 1 — ~ 
Ai - A 2 + A 1 (l + z) + A 2 (l + z) 2 to study the evolu- 
tions of DE EoS. The baryon acoustic oscillation (BAO) 
measurement from SDSS and the shift parameter de- 
termined from WMAP3 data combined with the new 
182 gold SN la data are used in our analysis. In sec- 
tion IV, we fit the simple one-parameter representation 
w{z) = wo exp[z/(l + z)]/ (1 + z) to the combined SN la, 
SDSS and WMAP3 data to obtain the geometry of the 
Universe. Note that the simple one-parameter model fits 
the observational data as well as the two-parameter mod- 
els do. In section V, we conclude the paper with some 
discussions. 
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II. RECONSTRUCTION OF THE 
DECELERATION PARAMETER 

The Hubble parameter H(t) — a /a and the decelera- 
tion parameter q(t) — ~d/(aH 2 ) are related by the fol- 
lowing equation, 



H(z) = H exp 



[l + q(u)]dln(l+u) 



(1) 



where the subscript means the current value of the vari- 
able. If a function of q{z) is given, then we can find the 
evolution of the Hubble parameter. For the flat ACDM 
model, q(z) = [fi m (l + zf - 2(1 - fi ro )]/2[fi ro (l + zf + 
1 — f2 m ] . In this section, we use the simple two-parameter 
function 



(2) 



to reconstruct the evolution of q(z). Note that go = 
1/2 + g 2j and dq/dz\ z=0 = qi — 2g 2 , so the parameter q 2 
gives the value of go- At early times, z ^> 1, q(z) — > 1/2. 
Substitute Eq. © into Eq. (JIJ, we get 



ff(«) =iJ (l + 2) 3/2 exp 



32 

2 



01 Z 



92 



2(1 + 



(3) 



Since the expression for the Hubble parameter is explicit, 
so we can think that we are actually modelling H{z) in- 
stead of q{z). 

The parameters q\ and q 2 in the model are determined 
by minimizing 



E 



(4) 



where the extinction-corrected distance modulus fi(z) = 
51og 10 [di(z)/Mpc]+25, <Ji is the total uncertainty in the 
SN la data, and the luminosity distance is 




FIG. 1: The evolution of q(z) = 1/2 + ( qi z + q 2 )/(l + zf by 
fitting it to the 182 gold SN la data. The solid line is drawn 
by using the best fit parameters. The shaded areas show the 
la, 2a and 3a errors. 



table HI we list the number N of SN la in a given rcdshift 
range for the 182 gold SN la data. The behavior was also 
found in 0,013 El BE, El in the fitting of the EoS 
of DE for a variety of models. This may suggest that the 
behavior of DE can be better constrained in the redshift 
range 0.1 < z < 0.6. The sweet spot can be estimated 
from the covariant matrix of errors, which is the inverse of 
the Fisher matrix in the linear approximation [HI, EH, El • 
The Fisher matrix is estimated to be Fu — 2.37, F\2 — 
F21 = 9.55, and F22 = 43.9. By choosing a\ = q\ and 
Q'2 = (Fi 2 /F 22 )qi +<Z2 = 0.2175gi + (72, the Fisher matrix 
becomes diagonal, ct\ and ct2 are uncorrelated, and the 
errors of ai and 0.2 are a 2 (ai) = -^22/(^11^22 — F? 2 ) and 
<7 2 (a2) = -F22 1 [37]- In terms of a\ and (X2, we get 



?(*) 



1 ai (z - 0.2175) +a 2 

2 + (1 + z) 2 ' 



(6) 



d L {z) = {l + z) 



dz' 

w~y 



(5) 



Fitting the model to the 182 gold SN la data, we get 
X 2 = 156.25, q-i = 1.47±^|| and q 2 = -1.46 ± 0.43, here 
the given error is the la error. For comparison, we also 
fit the ACDM model to the 182 gold SN la data and find 
that x 2 = 156.16, Q m = 0.48±g;^ and Q k = -0.44^°;^. 
So the simple two-parameter model of q(z) fits the SN 
la data as well as the ACDM model does. By using the 
best fitting results, we plot the evolution of q(z) in Fig. 
EJ From Fig. [2 we see that q{z) < for < z < 0.2 
at the 3cr confidence level. This result is consistent with 
previous analysis by using the 157 gold SN la data [l(|. It 
is also interesting to note that the stringent constraint on 
q(z) happens around the redshift z ~ 0.2. One may think 
perhaps there are more SN la data around this redshift. 
On the contrary, less SN la data is around z = 0.2. In 



Now the sweet spot can be estimated from the following 
equation 



2[ f 7 2 (a 2 )+a 2 (a 1 )(z-0.2175) 2 ] 
1 + z 



The sweet spot is estimated to be 



cr 2 (ai)(z- 0.2175). 

(7) 

t 0.2175 since 



a 2 (a 2 ) = F 2 . 



0. For a general model w(z) — w\ 



W2f(z) with arbitrary function f(z), the sweet spot is 
determined similarly from the equation f(z) — F12/ F\\. 



III. DARK ENERGY PARAMETRIZATION 

In this section, we use the observational data to re- 
construct the EoS of DE. For simplicity, we consider the 
spatially flat case, k = 0. In addition to the SN la data, 
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TABLE I: The distribution of the gold SN la data 



z 


< 0.1 


0.1-0.2 


0.2-0.3 


0.3-0.4 


0.4-0.5 


0.5-0.6 


0.6-0.7 


0.7-0.8 


0.8-0.9 


0.9-1.0 


> 1.0 


N 


36 


4 


5 


12 


31 


22 


16 


11 


17 


12 


16 



we also use the distance parameter 



A = 



0.35 



0.35 
£(0.35) 



0.35 



dz 

W) 



1/3 



(8) 



measured from the SDSS data to be i = 
0.469(0.95/0.98)-° 35 ± 0.017 0, II, and the shift 
parameter (3lj | 



K = 



1.70 ±0.03, 



(9) 



where E(z) = H(z)/H and z ls = 1089 ± 1. 

The first DE parametrization we consider is [l7| 

w(z) — W + 



l + z 



(10) 



The dimensionless DE density is 

n DE (z) = n DE0 (i + zf( 1+wo+w ^ 



exp[ 



-3w a z/(l + z)}. 

(11) 

This parametrization can be thought as the parametriza- 
tion of the DE density instead of w(z). By fitting this 
model to the observational data, we find that \ 2 — 



158.07, n m = 0.29 ± 0.04, w 



-1.07 



+0.33 
-0.28 



and w a 



0.85_^3 8 . Compared with previous fitting results [16l |. 
the current data makes a little improvement on the con- 
straint of w a . The evolution of w{z) is plotted in Fig. 
[2] and the contours of wq and w a are shown in Fig. [3] 
From Fig. [21 we see that at the 3a confidence level, 
w(z) < for z < 2, w(z) crosses the —1 barrier around 
z ~ 0.1, and the stringent constraint on w(z) happens 
around z ~ 0.3. From the Fisher matrix estimation, we 
get z/(l + z) — F12/F11 — 0.32, so the sweet spot is 
around z — 0.47. This estimation is quite different from 
what we get, and the main reason is that the distribution 
of w a is highly non-Gaussian. From Fig. [31 we see that 
the cosmological constant is more than la away from the 
best fit result. 

The second DE parametrization we consider is [l8j 



w(z) = Wq + 



w a z 



(12) 



(l + z) 2 ' 
The dimensionless DE density is 

n DE (z) = n DE0 {l + z) 3 ( 1+ ^) exp [3w a z 2 /2(l + zf] . 

(13) 

Again this parametrization can also be thought as the 
parametrization of the DE density. By fitting this model 
to the observational data, we find that x 2 = 157.11, 
Q m = 0.28t°;°t WQ = _i.37+o-B8 and Wa = 3.39l 3 ;^. 




FIG. 2: The evolution of w(z) by fitting the model w(z) = 
Wo + w a z/(l + z) to the observational data. The solid line 
is drawn by using the best fit parameters. The shaded areas 
show the la, 2a and 3<r errors. 




FIG. 3: The la, 2a and 3a contour plots of wo and w a for 
the model w(z) = wo + w a z / (1 + z) . The diamond denotes 
the point corresponding to the cosmological constant. 



These constraints are almost at the same level as previous 
results in [lj|. The evolution of w(z) is plotted in Fig. [4] 
and the contours of wq and w a are shown in Fig. [5l From 
Fig. [H we see that at the 3a confidence level, w(z) < for 
z < 0.7, w(z) crosses the —1 barrier around z ~ 0.15, and 
the stringent constraint on w(z) happens around z ~ 0.2. 
The sweet spot is estimated to be z — 0.2353 from the 
equation z/(l + zf = F 12 /F n = 0.1542. From Fig. [3 
we see that the ACDM model is almost 2a away from the 
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best fit result. 




Redshift z 



FIG. 4: The evolution of w(z) by fitting the model w(z) — 
Wo + w a z/(l + z) 2 to the observational data. The solid line 
is drawn by using the best fit parameters. The shaded areas 
show the la, 2a and 3a errors. 




FIG. 5: The la, 2a and 3a contour plots of wo and w a for 
the model w(z) = wo + w a z/(l + z) 2 . The diamond denotes 
the point corresponding to the cosmological constant. 



The last model we consider is [14| 

n DE (z) = A 1 (i+z)+A 2 (i+ z ) 2 +i-n r 

The EoS parameter w{z) is 

l + zAi + 2A 2 (l + z 



-A 1 -A 2 . (14) 



w{z) 



1. 



(15) 



The cosmological constant corresponds to A\ = A 2 = 0. 
By fitting this model to the observational data, we find 
that x 2 = 158.48, O m = 0.30 ± 0.04, A 1 = -0.48±^ 
and A 2 = 0. 25+^45. Tne evolution of w(z) is plotted 
in Fig. [5] and the contours of A\ and A 2 are shown in 



Fig. [7] From Fig. [6l we see that at the 3er confidence 
level, w(z) < for z < 1.1 and the stringent constraint 
on w(z) happens around z ~ 0.4. From Fig. [71 we see 
that the ACDM model is more than la away from the 
best fit result. Comparing the value of \ 2 °f t ne three 
models we considered, we find that the second model fits 
a little bit better than the other two models do. 




Redshift z 



FIG. 6: The evolution of w(z) by fitting the model £Ide{z) = 
1 - fi m - Ai - A 2 + A! (1 + z) + A 2 (1 + zf to the observational 
data. The solid line is drawn by using the best fit parameters. 
The shaded areas show the la, 2a and 3a errors. 




FIG. 7: The la, 2a and 3a contour plots of Ai and A 2 . The 
diamond denotes the point corresponding to the cosmological 
constant. 



IV. THE GEOMETRY OF THE UNIVERSE 

In this section, we use the observational data to find 
out the geometry of our universe. When k ^ 0, the 
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luminosity distance becomes 
dh{z) 



H y/W\ 



dz' 



E(z') 



(16) 



where sinn(y \k\x)/W \k\ = sin(x), x, sinh(a;) if k = 1, 0, 
— 1, the parameter A becomes 



0.35 



0.35 1 
£(0.35) |^ 



and the shift parameter becomes 



K 



ismn I \/\ttk, , 

^\ V 1 k, J E{z) 



ib dz 
o E(z) 



dz 



1/3 



(17) 



(18) 



To fit the observational data, we consider the one pa- 
rameter DE parametrization [la ] 



w(z) = - ' 



W r z/(l+z) 
l + Z 



(19) 



During both the early and future epoches, w{z) — > 0. 
The DE density is 



DE 



«dbo(1 + zf exp (3w e z/(1+z) - 3oj ) • (20) 



By fitting this model to the observational data, we find 
that x 2 = 158.85, TO = 0.30 ±0.04, Q k = -O.OOOT+q 1 ;^ 2 
and wo = — 0.93 jlp '^g. The data is also used to fit 
the ACDM model, the results are x 2 = 160.51, tt m = 
0.30 ± 0.03 and tt k = -0.02 ± 0.02. This model fits the 
observational data as well as the ACDM and the two- 
parameter models do. The contours of fl m and fl k are 
shown in Figs. [8]and[9l The errors on TO and are 
almost the same. Comparing with the results in [16| . we 
find that the new SN la data improves the constraint on 
0^ significantly. 




FIG. 8: The la, 2a and 3a contour plots of Q m and Qk for 
the parametrization w(z) — wo exp[z/(l + z)]/(l + z). 




FIG. 9: The la, 2a and 3a contour plots of Q m and Qk for 
the ACDM model. 



V. DISCUSSION 

By fitting the simple two-parameter representation of 
q(z) to the new 182 gold SN la data, we find strong ev- 
idence of acceleration in the recent past which is con- 
sistent with previous studies in 0, While the ev- 
idence of current acceleration is weak from fitting the 
simple piecewise constant acceleration model to the pre- 
vious gold SN la data [9J and fitting the simple two- 
parameter representation of q(z) to the 115 nearby Su- 
pernova Legacy Survey (SNLS) SN la data , we find 
strong evidence of current acceleration by using the gold 
SN la data. The strongest evidence of acceleration again 
happens around the redshift z ~ 0.2. The transition red- 
shift when the Universe underwent the transition from 
deceleration to acceleration is found to be z t — 0.361° ^ 
at the la level. 

The new SN la data, together with the BAO mea- 
surement from SDSS and the shift parameter determined 
from WMAP3 data, are used to fit three popular DE pa- 
rameterizations. When we are given the parameteriza- 
tions w{z) = WQ + w a z/(l + z) and w(z) = wq + w a z/{l + 
z) 2 , the explicit analytical expressions for the DE density 
can be derived. Alternatively, we can think that we are 
parameterizing the DE density Qoe{z) instead of w(z). 
The new observational data makes slightly improvement 
on the constraint of w a , while the ACDM model is still 
consistent with current observational data. Although 
high redshift SN la data provides robust constraint on 
the property of DE [32j], the stringent constraint on w(z) 
happens around z ~ 0.3. In other words, at the 3a con- 
fidence level, w(z) is best constrained around the red- 
shift z ~ 0.3. Surprisingly, we only have a few SN la 
with redshift around 0.3 in the current 182 gold SN la 
data. The same result holds for the DE parametrization 
fl DE (z) = l-n m -A 1 -A 2 + A 1 (l + z) + A 2 (l + z) 2 , al- 
though the redshift is now around z ~ 0.4. We think this 
result is quite generic for two-parameter parametriza- 
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tions. The result also suggests that more SN la data 
with the redshift z ~ 0.2 — 0.4 may be valuable to give 
strong constraint on w(z). The sweet spot around the 
redshift z <~ 0.2 — 0.4 may be argued from the Hubble 
law and the decreasing importance of DE [l^, [H, [HI : (1) 
At low redshift, the luminosity distance can be expressed 
as HocIl(z) = z + 4(1 — qo)z 2 . To the linear approxima- 
tion, it does not depend on the cosmological parameters, 
so the constraint on the property of DE at low redshift 
from SN la data is not strong; (2) At high redshift, the 
role of DE diminishes. Depending on the model, the ev- 
idence for w(z) < at high redshift is different. 

We also apply the one-parameter parametrization 
w(z) = w exp(z/(l + z))/(l + z) to study the geom- 
etry of the Universe. Although the SN la data alone 
does not provide valuable constraint on the geometry, the 
new combined data improves the constraint on signif- 



icantly. At the 3cr confidence level, we have |f2fc| < 0.05 
for the model w(z) — woexp(z/(l + z))/(l + z) and 
-0.064 < D, k < 0.028 for the ACDM model. It should 
be stressed that the effect of the heterogeneous nature of 
the gold SN la data on the systematics is also important 
and it may impose potential problem when combining 
with WMAP3 data [34], [3f| [36[ . The more ho mog eneous 
SNLS SN la data avoids this problem [34|, |35|, [36j . 
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